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The emission volume sizes of pions and Kaons, r 7T ± Jr ± and r K ± K ± , measured in 
the hadronic Z° decays via the Bose-Einstein Correlations (BEC), and the recent 
measurements of taa obtained by through the Pauli exclusion principle are used to 
study the r dependence on the hadron mass. A clear r 1T ± 7I ± > r K ± K ± > r^A 
hierarchy is observed which seems to disagree with the basic string (LUND) model 
expectation. An adequate description of r(m) is obtained via the Heisenberg uncer- 
tainty relations and also by Local Parton Hadron Duality approach using a general 
QCD potential. These lead to a relation of the type r(m) ~ Constant / y/m. The 

present lack of knowledge on the / O (980) decay rate to the K°K channel prohibits 
the use of the r K o K o in the r(m) analysis. The use of a generalised BEC and I- 

spin invariance, which predicts an BEC enhancement also in the K^K° and tt^tt 
systems, should in the future help to include the r K o K o in the r(m) analysis. 



1 Introduction 

For more than three decades the Bose Einstein Correlations (BEC) of identical 
bosons, mainly charged pions, were utilised to estimate the dimension of their emis- 
sion source. In nucleus- nucleus reactions a clear dependence of r X7r on the size of 
the nucleus is observed. In e + e~ annihilations leading to hadronic final states no 
clear evidence is seen for a change in r n7r , the emission volume size, as the ^/s ee 
increases from 3.1 GcV (the J/tp mass) to the Z° mass (see Fig. [j]). In particular if 
one remembers the large variety of experimental procedure adopted by the different 
experiments. Nevertheless, in order to be sure to isolate the effect of the hadron 
mass we concentrated our study to e + e~ annihilation at one energy only namely, at 
the Z° mass. At this energy the very large statistics of hadronic Z° decays of high 
average charged multiplicity, accumulated by the four LEP experiments, ALEPH, 
DELPHI, L3 and OPAL, affords an excellent opportunity to investigate whether 
the dimension r of two-particle emitter size does depend on their mass. 
In our analysis we rely on the 1-dimcnsional BEC studies where the emission source 



"Invited talk given by G. Alexander at the XXIX Int. Symp. on Multiparticle Dynamics, 9-13 
August 1999, Providence RJ, USA. (to be published in the proceedings of this conference) 
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Figure 1: A compilation of r 7I ± 1T ± measured in e+e annihilations. 



is taken to be a sphere with a Gaussian density distribution 
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'Q ). Here qi are the 4-momenta 



of the two identical bosons, A is the strength of the BEC effect and r is the emission 
volume size. 



2 Measurements of the source size r 

2.1 The r^i^i and r K ± K ± 

In Table 1 we present the weighted average value of the r v ± 7T ± measured in the 
Z° decays where the large systematic error reflects the spread of the individual r 
values due to the choices of the reference sample. In the same table are also listed 
the r K ± K ± derived by DELPHI from the BEC of charged Kaons and the preliminary 
value of OPALL4 . 



2 



2.2 The rjj-ojj-o and the Generalised BEC 

The three LEP experiments, OPAL, DELPHI and ALEPH, have also measured the 
BEC of the K%K% pairsl obtaining the r K o K o values listed in Table 1. The K° S K° S 

pairs have three sources where two of them, the K°K° and K° K , are identical 
di-bosons. As for the boson-antiboson pair K°K , no BEC effect is expected if one 
sums up all their possible decay pairs KgKg, K^Kg, KgKl and K\K\. However 
by selecting only the KgKg or the pairs BEC enhancement is expected with 

a matching decrease of the KgK^ pairs. Thus the r K o K o could have been used for 
our r(m) investigation was it not for the fact that the upper side of the width of 
the enigmatic / O (980) resonance is above the KK threshold. Recently it has been 
shown u that by invoking the generalised Bose-Einstein statistics and using I-spin 
invariance one can evaluate the minimum BEC contribution of KgKg pairs to the 
low mass enhancement. If the over-all hadronic final state of the Z° decay is in a 
pure 1=0, here denoted as ip a , then one has the relation: 

Ex_P&o - K s (p)K s (-p)X] = (l/2)-£x P[% - K^ffiK^-fiXl+Ex PVPo - 
K°K X — > Kg{p)Kg(—p)X] where X represents the hadrons accompanying the 
Kaon pair. From this relation follows that there will exist a BEC enhancement 
in the K S K S pairs the height of which will at least be half of that present in the 
charged identical Kaon pairs. Since this analysis has not been done so far, the 
current quoted r K o K o values cannot serve our analysis. Here is to note that the 
generalised Bose statistics predicts that BEC enhancement should also occur in the 
7r^~7r® (not studied so far) system in a similar strength to that seen in the 
pairs. 

The r values obtained from the tt ± tt ± and K ± K ± BEC analyses of DELPHI and 
the average values obtained from existing LEP experiments are shown in Fig 3. A 
first indication for the decrease of r with the increase of the hadron mass is seen 
but the data within the errors are still compatible with a constant r value, or even 
with a slight rise, as m increases. This situation calls for a BEC analysis of higher 
mass hadrons like e.g. the rf (958) meson. However with its very low production 
rate of about 0.14 per Z° decay, this avenue is closed to us. 

2.3 Estimate o/taa 

It was pointed out some time ago by G. Alexander and H.J. LipkinS that r values 
for pairs of baryons can be measured by observing the onset of the Pauli exclusion 
principle as the two identical baryons approach their threshold in their centre of 
mass system. The ratio £, defined as £= (S = l)/[(S = 1) + (S = 0)], measures the 
relative contribution of the S=l state to the two A system. Experimentally this 
can be measured through the angular distribution of the angle between the decay 
protons of the two identical hyperons defined as y* = cos(pip 2 ) . Here pi and 
P2 are the momenta of the protons coming from the two A's decay defined after 
two transformations: first the A's are transformed to their centre of mass system 
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Table 1: BEC results obtained at LEP1. (a) From spin composition analyses where no need exists 
for reference sample and Coulomb correction. 
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r (fm) 


Experiment 


7r±7r± 
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0.78 ±0.01 ±0.16 


LEP Average 


K±K± 


0.82 ±0.22 ±°-*l 


0.56 ±0.08 ±° ° 7 6 


OPAL (99) 




0.82 ±0.11 ±0.25 


0.48 ±0.04 ±0.07 


DELPHI(96) 


KgKg 


1.14 ±0.23 ±0.32 


0.76 ±0.10 ±0.11 
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0.96 ±0.21 ±0.40 


0.65 ±0.07 ±0.15 


ALEPH(94) 




0.61 ±0.16 ±0.16 


0.55 ±0.08 ±0.12 


DELPHI(96) 


AA 




0.11 ±0.02 ±0.01 


ALEPH(99) 


AA( a > 




0.19 ±£o7 ±0.02 


OPAL (96) 






0.11 t°;°3±o.oi 


DELPHI(98) 






0.17 ±0.13 ±0.04 


ALEPH(99) 



and in the second each proton is transformed to the centre of mass of its parent 
A. Using this method OPAL, DELPHI and ALEPH have estimated the taa value 
which are listed in Table 1. In order to determine the combined taa value of these 
three experiments we have considered the values of £(Q) shown in Fig. ||a. These 
values were fitted to the expression £(Q) = 0.75 [1 — e^ rAA< 3 '] with the result 
shown by the solid line in Fig. 2a corresponding to 

r AA = 0.15tg:82 fm. 

In Fig. 2b we show the distribution Ax 2 = X 2 — Xmin with Xmin = 2>9- 

The DELPHI values (circles) and the LEP averaged values (triangles) of r„± n ±, 

r K ± K ± an d r AA are plotted in Fig. ^ where a clear hierarchy is seen, namely 

rv±„.± > r K ± K ± > r AA . 

The error bars represent the statistical and systematic errors added in quadrature. 
The ALEPH collaboration has also attempted to measure taa by considering, sim- 
ilar to the BEC studies, the ratio between the frequency of data A pairs to that 
of Monte Carlo generated sample as a function of Q. The observed depletion of 
pairs at low Q was used to estimate r A \ = 0.11 ± 0.02 ± 0.01 fm (see Table 1). 
This method, which has a relative small statistical error in comparison to the spin 
composition measurement, rests however on the assumption that the depletion of 
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Figure 2: (a) The S = 1 fraction, e, of the AA(AA) pairs measured as a function of Q by the 
LEP collaborations. The solid line represents the fit results as described in the text, (b) The 



Ax 2 



dependence on taa obtained in the fit to the combined LEP data. 



events at low Q is indeed due to the Pauli exclusion principle and not for any other 
reason. This method, furthermore, has to rely on a reference sample. 



3 Comparison between r(m) and theoretical models 

The fact that r{m) decreases with the mass of the particle can be checked against the 
predictions of multi-particle production models of e + e~ annihilations. In particular 
this r(m) behaviour seems to contradict the expectation of the string (LUND) model 
where r(m) should increase with the particle mass. 

On the other hand the r(m) behaviour can be described in terms of the Heisenberg 
uncertainty relations Q i.e., ApAr — 2 p r = m v r = he and AEAt — 
2[p 2 /(2m)] At = h where m is the mass of the hadron, v and p its velocity and 
momentum and r is spatial separation between the two hadrons. In AE, possible 
contributions from a potential energy are neglected. From these equations one 
obtains 

r(m) = cVhAl/y/m~ ~ 0.243/ y/m(GeV) fm 

when considering only the kinetic energy, equating Ar = r and setting At to 10~ 24 
seconds representing the strong interaction time scale. As can be seen in Fig. || 
the thin solid line, which represents the expectation from the uncertainty relations, 
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Figure 3: The measured emitter radius r as a function of the hadron mass m compared to some 
theoretical predictions (see text). 

describes fairly well the data both in magnitude and in shape. To illustrate the 
sensitivity to our particular choice of At we also show by the dashed lines the 
r(m) expectations for At = 1.5 x 10~ 24 and 0.5 x 10~ 24 seconds. Rearranging the 
expression for r(m) it is amusing to find out that mr 2 (m) = Const. ~ 0.06 GeV 
fm 2 i.e., the identical hadron pairs seem to emerge with the same moment of inertia 
irrespective of their mass value. n 
Another approach to describe the r(m) dependence has been tried ouu. Namely, 
by using the virial theorem in the frame work of the Local Hadron Parton Duality 
hypothesis. In this context the general QCD potential given by 

4 a s fic 

V (r) = k r — 

w 3 r 

was used with the parameters set of k — 0.14 GeV 2 = 0.70 GeV/fm with a s — 
2n/9 \n(S+j/r), 5 — 2 and 7 = 1.87 GeV -1 = 0.37 fm as obtained from the hadron 
wave functions and decay constants. The result of this QCD based approach is 
shown in Fig. 3 by the solid thick line differs only very slightly from the uncertainty 
relations outcome. 
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4 Conclusion and summary 

The source dimension r is found to change from w 0.75 fm for charged pions to w 
0.15 fm for A's. This trend is opposite in direction to that expected in the basic 
string (LUND) model where dr(m)/dm > 0. It can however be described in terms 
of the Heisenberg uncertainty relations as well as in the frame work of the Local 
Parton Hadron Duality assumption. Both expect that r(m) ~ Constant/ y/rn. 
This relation leads to the conclusion that pairs of hadrons are produced with the 
same moment of inertia irrespective of their mass value, which presents a challenge 
for future attempts to formulate multi-particle production models. More precise r 
measurements for the source dimension will be beneficial for the comparison with 
models. Of particular interest will be the result from a BEC analysis of the 7777 
system which has a mass near the K-meson but it is strange-less and its production 
rate in the Z° hadronic decays is relatively high. At the higher mass end there is 
little chance to further explore the r(m) behaviour. The production rate of the 
7y'(958) is too low and so are also the rates for the hyperons above the A baryon. 
Presently one cannot utilise for the r(m) analysis the measured r K o K o values due 
to the fact that they may be influenced by the enigmatic / D (985) resonance decay 
into the K°K channel. In the frame work of the generalised BEC and I-spin 
invariance there is a possibility to separate the two sources of the Kg Kg system. 
This generalised BEC extension can be verified by observing BEC in the n ± 7T° pairs 
coming e.g. from hadronic Z° decays. 
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